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Abstract 

The probability density function (PDF) of flux R is computed in sys- 
tems with logarithmic non-linearity using a model non-linear dynam- 
ical equation. The PDF tails of the first moment flux are analytically 
predicted to be power law. These PDF tails are shown to be broader 
than a Gaussian distribution and are a manifestation of intermittency 
caused by short lived coherent structures (instantons). 
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I Introduction 



One of the remaining challenges in classical physics is to understand the na- 
ture of turbulence. In particular, gaining information of turbulent systems by 
experimental determination of probability density functions (PDFs) of flux 
is very common. While traditional mean field theory is built on Gaussian 
statistics (such as transport coefficients), avalanche- like events cause the de- 
viation of the PDFs from a Gaussian prediction. In particular, PDF tails due 
to rare events of large amplitude are often found to be substantially different 
from Gaussian although PDF centers tend to be Gaussian. For example the 
the scaling of PDF tails of fluxes in tokamaks is often found to be exponen- 
tial [1]- jl]. In this paper we present a statistical model of intermittency by 
computing PDF tails. Our key method is based on the idea that coherent 
structures such as vortices, streamers and blobs cause intermittency in the 
PDF tails. Significant flux can be mediated by coherent structures through 
the formation of avalanche like events of large amplitude, as indicated by 
recent numerical simulations and experiments [5]- [TB] . 

Specifically, the theoretical method used here to compute the PDF tails 
is a non-perturbative technique, the so-called the instanton method [16J- 
|27j . The instanton method has been adopted from quantum field theory 
and then modified to classical statistical physics for Burgers turbulence [16J- 
[TT] and in a model by Kraichnan [28]. Note that there is another method 
determining the PDF tails, namely to compute the Fokker-Planck equation 
for the PDF [27], [21]. 

In previous papers it has been shown using the instanton method that 
the PDF tails of momentum flux and heat flux in plasma turbulence are 
significantly enhanced over the Gaussian prediction [18]- [2T]. In particular, 
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a novel explanation for the exponential scaling of momentum flux R of the 
form ~ exp{— c(i?/i?o) 3//2 } found in experiments on CSDX at UCSD has 
been provided [3J. It was also shown that shear flows can significantly reduce 
the PDF tails of Reynolds stress and zonal flow formation [22]- [22]- I n 
a cubic non-linear model of shear flows the PDF tails were found to be 
~ exp{— c(i?/i?o) 4 } and thus fall off much faster than a Gaussian, which was 
confirmed by numerical simulations [24]. 

Current models fall short in predicting or interpreting experiments and 
simulations of an inherently strong intermittent nature. In some cases this 
involves models that have exponential or even logarithmic non-linear inter- 
action. For example exponential non-linearities may be found in models of 
current sheaths [29] in magnetically confined plasmas. On the other hand, 
logarithmic non-linearities are found in a variety of models, including a log- 
arithmic non-linear Schrodinger equation (LNLSE) introduced in Refs. [30J- 
[31]. The LNLSE is appealing from a mathematical point of view in that it 
supports solitary wave solutions (Gaussons) while still preserving many sim- 
ple features of linear equations. Specifically in the LNLSE, non-interacting 
sub-systems are still separable and there exist a lower energy bound. The 
LNLSE has also been studied in a stochastic reformulation of quantum me- 
chanics (QM) where to each quantum state there exists a stochastic process 
determined by a Langevin equation [32J. However, in QM the non-linear 
effects have been found in experiments to be very small with a |c| ~ lQ~ l5 eV 
( |c| is defined below in Eq. (30)) [33]- [35]. Nevertheless, this type of equa- 
tion is used in many other areas of applications such as nuclear physics, 
optics and geophysics; see e.g. [36]- [39]. In quantum field theory logarith- 
mic non-linearities naturally appear in super-symmetric field theories and 
inflation cosmology; a Klein-Gordon equation with a logarithmic potential 
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was studied in Refs. (ID]- [2] and have direct applications of the models in 
Ref. 02]- gl]. 

In principle there are roughly four different types of non-linear inter- 
action terms for fluctuations: exponential non-linearities; power law non- 
linearities; logarithmic non-linear terms; and cyclic non-linear terms. How- 
ever for weak fluctuations all these could easily be transformed into power 
law non-linearities by Taylor expansion which result in PDF tails of the form 
~ exp{— c(R/Rq) s } with s = (n + l)/m. Here, n and m are the order of 
the highest non-linear interaction term and moments for which the PDFs 
are computed, respectively [23]. For arbitrary fluctuations, in the case of 
exponential non-linear interaction the PDF tails was found to be described 
by the Gumbel distribution which represents a frequency distribution of the 
extreme values of the ensemble [26] . However it is questionable that a Taylor 
expansion of the non-linear interaction term would be valid for an instanton 
driven process and thus a more rigorous study is needed. 

The purpose of the present paper is to provide analytical expressions of 
the PDF tails of flux of first moment variables such as potential, density etc 
in systems with logarithmic non-linearities, i.e. PDF tails of the field itself. 
The predicted PDFs are shown to be power law for first moment flux, which 
are thus enhanced over Gaussian predictions. 

The paper is organized as follows. In Sec. II the model non-linear dy- 
namical equation is presented together with preliminaries of the path-integral 
formulation for the PDF tails. In Sec III the instanton solutions are calcu- 
lated and in Sec IV the PDF tails are presented. A discussion of the results 
and conclusion is given in Sec. V. 
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II Non-perturbative computation of the PDF 
tails for logarithmic non-linearity 

Here we present a statistical theory of PDF tails in systems with a logarithmic 
non-linearity. We follow the procedure developed earlier [IB]- [27] to compute 
PDF tails in systems governed by logarithmic non-linear interactions. To 
elucidate the particular features of the logarithmic non-linear interaction we 
start by computing the PDF tails of the reduced case, keeping only the 
logarithmic non-linear interaction (ln0) and later determine the PDF tails 
of the equation with a term 01n0 proposed in [30]- [31]. It has been shown 
earlier that PDF tails are rather insensitive to details of the dynamics and 
depend only on the dominant non-linear term [23J. Thus, we consider a non- 
linear equation in one spatial dimension describing the time evolution of the 
variable <fi for fluctuations governed by a logarithmic non-linear interaction 

^ + cln(0)-r / V 2 = /. (1) 

Here the term cln(</>) is the logarithmic non-linear interaction term, rj is a 
damping term and / is the forcing. For simplicity, the statistics of the forcing 
is assumed to be Gaussian with a short correlation time modeled by the delta 
function as 

(f(x,t)f(x',t')) = 6(t-t') K (x-x'), (2) 

where (/) = and k(x — x') = Koe^ x ~ x '^ L2 . The angular brackets denote the 
average over the statistics of the forcing /. The delta correlation in time was 
chosen for the simplicity of the analysis. In the case of a finite correlation 
time non-local integral equations in time are needed. 

We calculate the PDF tails of first moment flux M(<p) (density, potential) 
by using the instanton method. In general, M(0) is the m multiple product 



of 4> {mth moment) and can be denoted M(0) = P m (0). In particular, for 
the first moment P\{(j>) = ai Pi where a is a constant. The probability density 
function for the potential 4> can be defined as 

P(R) = (5{M(<f>(x = x )) -R)) = j d\e lXR h, (3) 

where 

/ A = (expHAM(0(x = x o )))>. (4) 

The integrand can then be rewritten in the form of a path-integral as 

I x = J V(f)V0e- s \ (5) 

Here DtpVcj) is the path-integral measure. The assumption of Gaussian statis- 
tics of the forcing is used to formally express the PDF tails in terms of a path 
integral [27] . 



Ill Instanton (saddle-point) solutions for a 
system with logarithmic non-linearity 

The path integral in Eq. (5) will be computed by the saddle point method 
to evaluate the effective action S\ by requiring that the variational first 
derivatives vanish for the optimum path (among all possible paths). The 
saddle-point solution of the dynamical variable 4>(x, t) of the form 4>(x, t) = 
F(t)v(x) is called an instanton if F(t) — at t — — oo and F(t) ^ at 
t = 0. The optimum path is associated with the creation of a short lived 
coherent structure (instanton). Thus, the bursty event can be associated 
with the creation of a coherent structure. Note that, the function v(x) here 
represents the spatial form of the coherent structure. By using the new 
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variable = v(l + F), we express the action in Eq. (5) by 
S x = -ij dtcxh (P + c 2 ln(l + F) + r]c 5 (F + 1)) 

+ i\ Cl c 3 J dt(l + F)5(t). (6) 

Here, F\ is the time dependent part of conjugate variable <fi = vF\ with 
coherent structure v. Note that, the sub-leading order (in A) term vlnvFi 
has been neglected. The coefficients ci,C2, C3, C4 and C5 are defined as 

ci = J dxv(x)v(x), (7) 

cic 2 = J dxcv(x), (8) 

cxc 3 = J dxv(x), (9) 

C1C4 = J dxdyv(x)v(y)K,(x — y), (10) 

Cl c 5 = / rfx(V 2 0)0. (11) 



The term rj(F + 1) is due to dissipation. To find the instanton solutions we 
compute the first variational derivatives of S\ and take them to be zero 
S S 

~ = -tc 1 (F + c 2 \n(l + F)+r ] c 5 (l + F))+c 1 c 4 F 1 = 0, (12) 

^ = -ici ^1 + c 2T ^ + ^ 5 F^ + iAc lC3 <$(f) = 0. (13) 

The equation of motion (EQM) for t < is derived by differentiating Eq. 
(12) and substituting the result in Eq. (13) 

£(*+*ttf + = (14) 

which leads to 

-F + c t^^P + C 5 V(1 + F) + c 2 c 5 77 ln(l + F) + c 2 c 5V = 0. (15) 
1 + r 

We assume that dissipation is small but finite and drop the last three terms 
in Eq. (15). This is valid if T} 2 \ and r\ In F are sufficiently small. We integrate 



Eq. (15) by using u = dF/dt and F = udu/dF and use separation of variables 



to find 



dF_ 
~dt 



±c 2 m(l + F). 



(16) 



The instanton function (F) now requires the evaluation of the following log- 
arithmic integral 



Here, F = lim e __>.o F(t = e_) is the boundary condition at t — and 
F t = F(t) for t > ^ where too is some large but finite negative time which 
is bounded due to dissipation. That is, the effect of damping cannot be ne- 
glected for t < (note t, too < 0). The resulting logarithmic integral would 
yield an instanton only when a finite value of dissipation is incorporated. In 
practice, it is important to note that to compute PDF tails from the action 
given in Eq. (6) we need only know the value of the first time derivative of 
the instanton (F) and its boundary condition at t — (F ). 

Upon integrating Eq. (13) over (— e, e), the boundary condition Fi(e) = 
gives 




(17) 



F 1 (-e) = c 3 A. 



(18) 



We then use Eq. (12) at t — to obtain 



ln(Fo) 



2C 3 C 4 



A. 



(19) 



2c 2 
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2ci 

c 4 
2cic 2 



IV PDF tails for the logarithmic non-linear 
interaction 

Here we compute the A dependence of the action S\ ~ h(X), which determines 
A- integral in Eq. (3). We first evaluate the action in the large A limit 

S x = -i f dt2c 1 c 2 F 1 \n(l + F) 
+ i Cl c 4 J dtFl 

+ iAcic 3 J dtF8(t) (20) 
J dt(\n(l + F)f + i\ Cl c 3 J dtF5(t) (21) 

(F(0) ln(F(0)) - F(0)) + iX Cl c 3 F(0). (22) 
c 4 

The action can only be estimated asymptotically for large A since there is a 
small constant contribution at too. The action S\ can be expressed as 

S x w -e- iaX , (23) 

where the parameter a is 

The tail of the PDF is found by computing the integral in Eq. (3) 

P(R) ~ y ' d\e lXR+e ~ ia \ (25) 

In order to use the saddle point method on the A integral, we let 

/(A) = iRX + e~ ia \ (26) 

and find A , which gives the extreme of the function / (i.e f'(Xo) = 0) as 
iR = iae~ iaXo . (27) 
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By taking the natural logarithm on both sides, we conclude that the resulting 
saddle point is 

Xo^ia-Hnia^R). (28) 

Using this into Eq. (25) gives us 

P(R) ~ e -# ln #+ e ' n# 

= e-^f-D. (29) 

The resulting PDF is thus a power law to leading order when R is large. 

Next we will consider the 01n0 non-linear interaction LNLSE and we 
will show that the PDF tails of first moment flux are close to a Gaussian 
distribution. In addition, we will calculate the PDF tails of the 2nd moment 
(\4>\ 2 ) flux and we will show that these tails are similar to the tails of first 
moment flux of the ln0 model (Eq. 29). Following the method used above 
we find the relevant dynamical equation to be 

</> + c</>ln|</>| 2 = /. (30) 

Here / is defined in the same way as above. The action can be written as 
(neglecting subleading order terms) 

S x « -if dtdF^F + 2cF In F) + J dtF? 

+ iXc lC3 j dtF k 5(t). (31) 

Here F is assumed to be positive and k — 1 the for first moment and k — 2 
for the second moment flux. We first consider the case k — 1. The first varia- 
tional derivatives of the action S\ give dynamical equations for the instanton 
function F and its conjugate F 1 that we use to find a relation between the 
function F and the time derivative F as 

F = 2cF In F, (32) 
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and an equation for the value of F at t = 



F(0)lnF(0) 



-«c 3 c 4 



A. 



(33) 



4c 



The action can now be computed in the large A limit to yield 




(34) 



16c 



The PDF tails are now found by integration using the saddle-point method 
with the saddle point A = —iR/(2k\) 



Here k\ = CiC3C 4 /(16c) and k 2 = l/(cic 3 ). We note that these PDF tails 
are very close to the tails of a Gaussian distribution. In the case of the 2nd 
moment flux |0| 2 (k = 2) we find 



with k 3 = — + ^ and k 4 = The PDF tails of 2nd moment flux in 

u C3C4 C4 ^ C4 

LNLSE are significantly enhanced compared to a Gaussian distribution. The 
PDF drom the 2nd moment are similar in structure to the PDF tails of the 
first moment in the In model due to the logarithmic non-linearity. 

V Discussion and conclusions 

In this paper we have presented the first computation of PDF tails of flux 
in dynamical system with logarithmic non-linear terms (ln0). Specifically, 
the PDF tails of the first moment flux was shown to be a power law. This 
result is rather insensitive to the precise form of dynamical equation and is 
only dependent on the dominant non-linear term [23] . The PDF tails deviate 
significantly from Gaussian distributions signifying that rare events with high 
amplitude influence the total flux. 




(35) 



P(R) 



-k 3 RlnR+k 3 R 



(36) 
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Figure 1: The PDF tails of flux in a system with logarithmic non-linearity 
and a Gaussian distribution are shown. The values of the parameter a are 
a = 1.0 (blue, dashed line); a = 0.5 (black, dotted line); a = 2 (green, solid 
line); and a Gaussian with a coefficient of unity (red, dashed-dotted line). 

We now show explicitly that systems with logarithmic non-linearities have 
fat PDF tails by plotting [Eq. (29)] and comparing the results to a Gaussian 
distribution. The shown distributions are normalized as one-sided distribu- 
tions. In figure 1, the PDF tails of flux with a logarithmic non-linearity and 
a Gaussian distribution are compared. In the figure the values of the param- 
eter a are a = 1.0 (blue, dashed line); a = 0.5 (black, dotted line); a = 2 
(green, solid line); and a Gaussian with a coefficient of unity (red, dashed- 
dotted line). For logarithmic non-linearity a stretched tail (a power law) is 
evident. This suggests that the resulting PDF tails in systems with logarith- 
mic non-linearities are stretched compared with Gaussian distributions and 
that rare bursty events may be predominant in transport processes. Note 
that the coefficient a is dependent on the coherent structure and the forcing 
through the constants C2, C3 and C4. 
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Moreover the usefulness of these PDF tails is limited by the spatial struc- 
ture of the specific solutions to the non-linear equations. The spatial struc- 
ture determines the magnitude and sign of the coefficient a. 

In summary, in the present paper we have investigated a statistical the- 
ory of turbulence and intermittency due to coherent structures in dynamical 
systems with logarithmic non-linearity. The use of coherent structures is 
motivated by various experimental results where bursty events cause a sig- 
nificant transport which are linked to coherent structures. We have computed 
the PDF tails of the 1st moment (e.g. density, potential) and 2nd moment 
(|0| 2 ) caused by intermittent coherent structures using the non-perturbative 
instanton method. The obtained power law PDF tails of first moment flux of 
ln0 model (Eq. (29)) are shown to be significantly enhanced in comparison 
with a Gaussian distribution. For the LNLSE model with 01n0 interaction 
we found PDF tails of the first moment flux close to Gaussian, however the 
2nd moment flux were similar to the PDF tails of 1st moment flux of the 
hvcj) model. The precise form of the PDF tails in non-linear Klein-Gordon 
system is however beyond the scope of the present paper and will be studied 
in future publications. Note that, although Gaussian forcing is used, the 
non-linearity in the system can give PDF tails that are far from Gaussian. 
Furthermore, the PDF tails are rather insensitive to the precise details of 
the dynamical equation and depend only on the dominant non-linear term. 
This result could provide an important piece of information for interpreting 
experimental results; the PDF tails derived in here and in previous papers 
may directly be compared to experiments. 
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